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Abstract 

Ionization and excitation of water molecules in intense laser pulses is studied theoretically by solving the three- 
dimensional time-dependent electronic Schrodinger equation within the single-active-electron approximation. The pos- 
sibility to image orbital densities by measurement of the orientation-dependent ionization of H2O in few-cycle, 800 nm 
linear-polarized laser pulses is investigated. While the highest-occupied molecular orbital 1 bi is found to dominate the 
overall ionization behavior, contributions from the energetically lower lying 3 ai orbital dominate the ionization yield in 
the nodal plane of the 1 61 orbital. The ratio of the ionization yields of the two orbitals depends on the intensity. Further- 
more, even for laser pulses as long as 8 cycles the orientation-dependent ion yield depends on the carrier-envelope phase. 
In the interpretation of the orientation-dependent ionization as an imaging tool these effects have to be considered. 

Keywords: Water, H2O, Multiphoton ionization, Density-functional theory. Orientation dependence. Orbital imaging, 
Ultrashort laser pulses. Time-dependent Schrodinger equation 
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1. Introduction 

Recent efforts to achieve time-resolved imaging of 
chemical reactions or other structural changes with sub- 
femtosecond resolution possess a strong interdisciplinary 
character, since the goal is of interest for physics, chem- 
istry, and biology. Pioneering experiments with near- 
infrared ultrashort laser pulses like the ones in [U [5] in 
which structural information about the chemically most 
relevant valence electrons was extracted from the emitted 
high-harmonic radiation or electron spectra have been of 
great interest, since these approaches should intrinsically 
have the potential for providing also the required time res- 
olution. However, in contrast to the initial assumptions, 
recent experiments have demonstrated that the molecular 
strong-field response depends, at least for some molecules, 
on more than one orbital |3]-E]- While such multi-orbital 
effects clearly complicate simple imaging schemes, they 
can also be the source for even richer information that 
can be gained from such experiments. An example is the 
electron-hole dynamics in the laser-generated ion that may 
be observed by analyzing the high- harmonic radiation [S]. 
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Though very exciting by itself, this appears to make di- 
rect imaging of the valence electrons and their field-free 
dynamics during, e.g., a chemical reaction more compli- 
cated. 



The proposed imaging schemes based on linear-polarized 
ultrashort near-infrared laser pulses may roughly be di- 
vided into two categories, rescattering-based schemes and 
direct imaging. The first category is based on the cele- 
brated three-step model of strong-field physics in which 
(1) an electronic wavepacket leaves the molecule around 
the local maxima of the electric field by tunneling ion- 
ization, (2) this wavepacket is accelerated in the laser field 
and reverses its direction as the field direction changes, and 
(3) the electronic wavepacket may recollide with its par- 
ent ion. As a consequence of this recoUision, the electronic 
wavepacket may partly scatter elastically (diffraction) or 
inelastically (leading to excitation or further fragmenta- 
tion), or recombine by the emission of high-harmonic ra- 
diation. Clearly, all these processes should depend on the 
structure of the molecular ion and thus have the potential 
for revealing structural information. This includes both 
electronic structure as well as nuclear geometry. Corre- 
sponding reviews may be found in [3 [S] . 
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2. Direct imaging using orientation-dependent 
ionization 

While the rescattering-based imaging schemes are evi- 
dently based on the third step of the three-step model and 
consequently require that this model is at least a good 
approximation, the direct imaging schemes are based on 
the first, the ionization step. Measuring the ion yield as a 
function of the time-delay between two very short pulses 
has been demonstrated theoretically [S] and experimen- 
tally [lOj to induce and image vibrational motion in the 
electronic ground state of neutral molecules. Simplified 
models like the strong-field approximation of the time- 
dependent Schrodinger equation suggest also the possibil- 
ity to image the nuclear geometry of a molecule by com- 
paring the energy-resolved electron spectra for different 
molecular orientations, but were immediately shown to fail 
using a full-dimensional ab initio treatment for H2 
An alternative, rather straightforward approach to time- 
resolved direct imaging with intense near-infrared laser 
pulses is based on recording the dependence of the total 
ionization probability as a function of the relative orien- 
tation of the polarization axis of a linear polarized laser 
pulse with respect to the molecule. Since this corresponds 
to the fully integrated electron (or ion) spectrum, it is by 
far the most intense response signal, many orders of mag- 
nitude larger than recollision signals. This makes it exper- 
imentally very attractive. In [12 it was demonstrated that 
the alignment-dependent ionization probability reflects the 
shape of the highest-occupied molecular orbital (HOMO) 
for both N2 and O2. However, the therein also considered 
example of CO2 turned out to be not as straightforward. 
A possible reason could be a coherent core trapping of the 
valence electron [13^ which explains that only at lower in- 
tensity the imaging appears to work for CO2 , in agreement 
with an experiment performed with lower laser intensity 

El- 

Finally, it should be reminded that time-resolution is 
achieved in direct imaging schemes using pump-probe ar- 
rangements. Typically, the dynamics has anyhow to be 
triggered in some way. This may be achieved by two iden- 
tical, time shifted laser pulses as in [51 [TUj or by synchro- 
nizing another pump laser to the probe pulse. On the 
other hand, especially the high-harmonic based recollision 
imaging schemes possess an intrinsic time resolution that 
is, however, limited in total duration to about half a laser 
cycle and thus to about 1.6 fs for 800 nm laser light [T3]. To 
follow dynamics on longer time scales, lasers with longer 
wavelengths have to be used, or a pump-probe scheme has 
to be applied as in the direct imaging schemes. Further- 
more, the intrinsic time resolution can only be used for 
processes that are induced in the ionizing first step and 
thus not really adequate for the general purpose of imag- 
ing chemical reactions. 

In order to image chemical reactions, it is, of course, im- 
portant that concepts tested so far for linear molecules are 
also applicable to more general molecules. We have thus 



decided to investigate the direct imaging proposal based 
on the recording of the orientation-dependent ionization 
yield for a simple but very important non-linear molecule: 
water. Despite its natural importance due to its abun- 
dance in nature, there are some additional arguments for 
choosing this molecule in the present context. The three 
energetically highest-lying orbitals of H2O are known to 
very roughly correspond to the three p orbitals of the oxy- 
gen atom that are, however, differently modified by the OH 
bonds. All three orbitals are thus structurally relatively 
simple and a good test candidate for imaging schemes. For 
example, the energetically highest-lying occupied orbital, 
the 1 bi HOMO, is almost identical to an oxygen p orbital. 
In contrast to true atomic p orbitals, the two hydrogen 
atoms in water lead, however, to a break of the isotropic 
symmetry and a deformation, especially of the so called 
3ai HOMO-1. Furthermore, due to the symmetry reduc- 
tion the three oxygen p orbitals loose their degeneracy and 
thus become (also energetically) distinguishable. As a con- 
sequence, it is possible to test the sensitivity of an imaging 
scheme. Note, it has to be mentioned that the field-free 
alignment of a dipolar molecule like H2O is experimen- 
tally a challenging task. However, orientation-resolved in- 
formation for one-photon ionization was obtained in [16] , 
but from a multiple coincidence measurement that allowed 
to reconstruct the orientation at the incident of ionization. 
Other approaches that avoid pre-alignment of the molecule 
are based on the use of circular polarized laser pulses as 
was applied to H2 in 17 . Clearly, in this case the theo- 
retical description should consider circular and not linear 
polarization. Due to the increased numerical demands, 
this is, however, so far seldom the case (see, e. g., the the- 
oretical analysis of the just mentioned H2 experiment in 
[HITTIIIS]). 

3. Previous studies of H2O in intense laser fields 

Ionization and harmonic generation of water molecules 
in intense fields have recently stirred some interest and 
were investigated both theoretically and experimentally. 
The case of relatively low intensities with supposedly neg- 
ligible ionization was studied using time-dependent con- 
figuration interaction (restricted to single and double ex- 
citations) in [19 where polarizabilities and low harmonics 
were reported. Ionization of the HOMO of H2O in electric 
half-cycle pulses was studied in [20^^ within a hydrogenic- 
orbital approximation and thus an atomic model. Adopt- 
ing a numerical implementation of time-dependent den- 
sity functional theory (TD-DFT) based on Voronoi cells 
the ionization of water molecules in intense laser fields 
was studied in j2T) . In that implementation only next- 
neighbor-cell interactions were considered and thus ap- 
proximate TD-DFT results were obtained. This approx- 
imate TD-DFT approach is in the following named TD- 
VFD as in [5T] . Considering two 2-dimensional cuts it was 
found that within the TDVFD approach the overall ion- 
ization is dominated by ionization from the HOMO, but 
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at specific directions the HOMO-1 contribution can locally 
exceed the one from the HOMO. In [52] the ionization from 
the HOMO of water was calculated within the molecular 
Ammosov-Delone-Krainov (MO-ADK) approximation. At 
least along the considered two-dimensional cut the depen- 
dence of the ionization rate on the relative orientation be- 
tween laser field and molecular plane as obtained by the 
MO-ADK model was shown to agree qualitatively, though 
not quantitatively, with the TDVFD results in |21) . 

The high-harmonic radiation of water molecules was re- 
cently investigated experimentally [23j and theoretically 
within a variant of the strong-field approximation in which 
the recombination dipole matrix element is related to 
single-photon ionization cross-sections ^24j. Already ear- 
lier, a theoretical investigation within the strong-field ap- 
proximation was reported in '25]. In accordance with the 
in ^21^ predicted overall dominance of the HOMO ioniza- 
tion (first step of the three-step model), the theoretical 
investigation of isotope effects on the high-harmonic gen- 
eration of water was restricted to the HOMO in j25j . It was 
concluded that practically no isotope effect is present and 
thus no nuclear motion is visible from the high-harmonic 
radiation, if the ratio of the high-harmonic spectra from 
H2O and D2O is analyzed. The reason is the absence of 
(substantial) excitation of vibrational motion, if ioniza- 
tion from the water HOMO is considered. This is due 
to the fact that the potential surfaces of the neutral wa- 
ter molecule and the formed ion (in its electronic ground 
state) are very similar. The time-resolved imaging of nu- 
clear motion using high-harmonics requires on the other 
hand vibrational excitation and thus different potential 
surfaces of the neutral and the ion ^ i26n28j . Neverthe- 
less, a recent mixed experimental and theoretical study [5] 
has shown a pronounced isotope effect in water harmon- 
ics. It is explained by high-harmonic generation from the 
HOMO-1 which, in contrast to the one from the HOMO, 
is accompanied by substantial vibrational excitation. As 
it becomes clear from the present study, the good visi- 
bility of the high-harmonic radiation from the HOMO-1 
(relative to the one from the HOMO) is, within the ter- 
minology of the three-step model, also due to the second 
and third step, while ionization (in an isotropic ensemble) 
occurs dominantly from the HOMO, at least for not too 
high intensities. 

The main goal of the present work is, however, 
to provide quantitative predictions of the full three- 
dimensionally resolved dependence of the ionization yield 
of a water molecule on the relative orientation of the 
molecule with respect to the laser field axis. Based 
on these results the prospects of direct (time-resolved) 
imaging of orbital densities by recording the orientation 
or alignment dependent ionization probabilities are dis- 
cussed. 




Figure 1: (Color online) Sketch (not to scale) of the spheres defining 
the local basis sets for H2O. The origin is chosen to coincide with 
the oxygen atom (blue), while the two H atoms (red) and thus the 
molecule lie in the xz plane. Two small spheres are centered around 
the hydrogen atoms, while one large sphere (with radius t^,^^ de- 
fines the main basis and discretization volume. The bold green arrow 
along the positive 2;-axis indicates the direction of the permanent 
dipole moment of water. 



4. Method 

The full ah initio treatment of the time-dependent re- 
sponse of molecules exposed to intense near-infrared laser 
fields is a non-trivial task. In fact, only very recently it 
has become possible to treat the two electrons of the sim- 
plest stable molecule, H2, exposed to intense ultrashort 
laser pulses (with 800 nm wavelength) for arbitrary ori- 
entations in full dimensionality [llj . Therefore, we have 
implemented a method that solves the time-dependent 
Schrodinger equation for in principle arbitrary molecules 
within the single-active-electron (SAE) approximation. 
The method was introduced in 29. where also the validity 
of the SAE approximation was demonstrated for H2 . The 
extension to linear molecules with more than two electrons 
(N2, O2, and CO2) followed in [13]. As it turned out, 
the standard frozen-core approximation is neither gauge 
invariant nor does it preserve the Pauli principle during 
the pulse, but it is possible to formulate a single-Slater- 
determinant approximation that does |30j . In the present 
work, almost all results were obtained with the latter ap- 
proximation (named method B in jl3j). Since details of 
the approach can be found elsewhere [131 [Ml 130], it 
is only briefly described with an emphasis on the present 
implementation for general, non-linear molecules. 

From the theoretical viewpoint the description of the 
ionization of non-linear molecules is challenging for two 
reasons. First, a proper description of the molecule has to 
account for its full multi-center character. Second, non- 
linear molecules belong to symmetry groups with no infi- 
nite symmetry element. Therefore, the use of symmetry 
does not lead to a dramatic reduction of the numerical 
efforts, since the Hamiltonian consists of a finite number 
of blocks and is by far not as sparse as for non-Abelian 
symmetries. The degree of sparseness determines on the 
other hand the efficiency of the time propagation. In the 
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18.69 
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30.50 


32.20 




535.16 


539.70 



Table 1: Calculated vertical ionization potentials {Ip) for the occu- 
pied orbitals of water compared to experimental values. 



case of H2O with symmetry there are only four ir- 
reducible representations: Ai, A2, Bi, and B2- In view 
of this, compared to linear molecules, small reduction of 
the computational efforts it is evident that the present 
study demonstrates the ability of our recently developed 
approach [T31 [51] to treat arbitrary molecules, because 
for the worst case, molecules with no symmetry (Ci), the 
number of irreducible representations is only reduced from 
four to one. 

In the first step of the calculation field-free Kohn-Sham 
orbitals are obtained using the LB94 exchange-correlation 
potential which has the correct asymptotic behavior |31| . 
The Kohn-Sham orbitals are described by a multi-centered 
i?-spline basis [32|. This basis is composed of a set of 
(typically atom-centered) spheres in which a local ba- 
sis is defined as the product of spherical harmonics and 
a radial part expressed in B splines. These spheres 
must (presently) not overlap. An arbitrary number of 
such spheres centered at any position of 3D space is al- 
lowed. Molecular symmetry is accounted for by generating 
symmetry- adapted basis sets. A large central sphere that 
is usually positioned at the charge center of the molecule 
is added that defines an additional large i?-spline basis. 
It overlaps with all atomic-centered spheres and provides 
an improved description of the chemical bonding and, es- 
pecially, the molecular electronic continuum. The latter 
is obtained in a discretized form, the density of the dis- 
cretized continuum being determined by the size of the 
central sphere. The concept is sketched for the present ex- 
ample of II2O in Fig. [1] Since the charge center lies very 
close to the oxygen nucleus, no local basis is defined for the 
latter, but the center of the large sphere coincides with the 
oxygen nucleus and not with the charge center. It should 
be emphasized that the adopted multi-center basis is well 
suited to treat the atomic Coulomb singularities and the 
atomic cores of a molecule. For example, in the frozen- 
core SAE approach in 33J a fitted core potential had to 
be smoothened (by hand) in order to prevent numerical 
instabilities arising from the Coulomb singularities. Fur- 
thermore, much higher angular momenta are required for 
describing a multi-center molecule in a single-center basis. 
Also the Cartesian coordinates used in [M] are not opti- 
mal, as the slow quantitative convergence of that approach 
demonstrates. 



In the present calculations, the radius of the large sphere 
is r'^^i^^ = 159.5 a. u. and the maximum number of angu- 
lar momenta is Zmax = 7. The experimentally determined 
equihbrium geometry (i?oH = 0.958i, <j:HOH = 104.5° 
[36]) is used and the water molecule lies in the xz plane. 
The calculated eigenenergies of the KS orbitals occupied in 
the ground state are given in Tab. [TJand compared to ex- 
perimental values taking the corresponding considerations 
in [3 7) into account. 

The time- dependent wavefunction is expressed as a sin- 
gle Slater determinant; the initial state being the ground- 
state Kohn-Sham Slater determinant. The one-electron 
wavefunctions that form the time-dependent Slater deter- 
minant are expanded in the basis of the occupied and 
virtual Kohn-Sham orbitals of the field-free Kohn-Sham 
molecular Hamiltonian. Thus the Kohn-Sham orbitals 
form a convenient basis, but in principle any other com- 
plete single-particle basis could be used, like, e.g., the 
Hartree-Fock basis additionally used in [29]. It should be 
emphasized that the present approach is not TD-DFT. 
In contrast, the orbitals are propagated in the (frozen) 
ground-state Hamiltonian plus the external field. As is 
shown in |30j , the time-propagation of the Slater determi- 
nant can be reduced to an independent time propagation of 
the initially occupied orbitals. The ionization yield is de- 
fined as the population of all field-free continuum orbitals 
after the pulse. The total population is normalized to 1, 
an ionization yield of, e.g., 0.1 corresponds to 10% ion- 
ization. As excitation we define the population of all field- 
free, initially unoccupied (virtual), but non-continuum or- 
bitals after the pulse. The inclusion of orbitals with an 
energy lower than 20 a.u. led to about 6,000 orbitals that 
were used in the calculations where the laser field is parallel 
to the x, y, or z axis due to symmetry and dipole-selection 
rules, while about 12,000 orbitals were used for arbitrary 
orientations. The laser pulses were modeled using a cos^ 
envelope for the vector potential, 

A{t) = ^0 cos2 (^^^ sin {ujt + 0) , (1) 

with pulse duration T, laser frequency tu, and carrier en- 
velope phase 0. The ionization yields are multiplied with 
a factor 2, since all orbitals of II2O are doubly occupied 
(and non-degenerate). This procedure was shown to be 
appropriate for H2, except for very high ionization yields 
(above about 10 to 20 %) where a factor 1 should be used 
|29j . The calculations were performed in length gauge and 
within the intervall ip G [0,7r/2] and e [0,7r], where ip 
and d represent the usual spherical coordinates. The full 
three-dimensional representations are obtained by symme- 
try arguments. 

5. Results 

The shapes (iso-surfaces of identical position-space den- 
sity) of the three energetically highest lying occupied 
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Figure 2: (Color online) Three dimensional plots of the normalized ionization (center column) and excitation (right column) yields of the 
HOMO (top), HOMO-1 (middle) and HOMO-2 (bottom). In each plot the silhouette (maximum circumference) is projected onto the yz 
(blue), xz (green) and xy plane (red) for a better visualization. The pulse parameters are a 800 nm, 8 cycle cos^ pulse with peak intensity 
2 X lO'^^ W/cm^. The shape of the individual orbitals (in real space) is given in the left column for comparison. 
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Figure 3: (Color online) Comparison of the cuts of the ionization 
yield of the HOMO-2 in xy and xz plane. The color coding agrees 
with the one in Fig. |2] 



molecular orbitals of H2O are visualized in the left panel of 
Fig. [2I Ordered by increasing electron binding energy, the 
orbitals are denoted as HOMO (1 61), HOMO-1 (3ai), and 
HOMO-2 (162)- The response of H2O to an 8-cycle cos^ 
pulse (see Eq. [T]) with a central wavelength of 800 nm, the 
laser peak intensity 2 x 10^'^ W/cm^, and = is also 
shown in Fig. |2] The ionization yields (middle column) 
and excitation yields (right column) are normalized with 
respect to the maxima, i. e. the maximum yield in each plot 
is one. In agreement with the shapes of the orbitals the 
ionization yields show approximately p-orbital like struc- 
tures that point along the directions of maximal extension 
of the orbitals. This is most evident for the HOMO that 
also itself is closest to an undistorted p orbital of an oxygen 
atom (the Pj^ in our chosen coordinate system). 

The HOMO-2 is also rather similar to an oxygen p or- 
bital (the Pa;), but the two lobes are distorted in the di- 
rection of the hydrogen atoms (within the xz plane) . The 
ionization yield of the HOMO-2 is not rotationally sym- 
metric around the y axis, but shows a similar distortion in 
the xz plane as the orbital itself. This can be seen even 
more clearly in the corresponding cuts shown in Fig. |3] in 
which the absolute values of the ionization yield as a func- 
tion of the orientation of the molecule with respect to the 
field axis is shown. However, the distortion is opposite to 
the one of the orbital, since it points toward the negative z 
axis, while the hydrogen atoms point toward positive z. In 
fact, one would not expect any asymmetry of this kind in 
in the strong-field response to a linear polarized laser field, 
since the electric field vector points periodically "up" and 
"down" . Therefore, the total ionization probability should 
be inversion symmetric. The found asymmetry is a conse- 
quence of the shortness of the adopted laser pulse, though 
8 cycles is not ultrashort. As a check, a calculation with 
a different carrier envelope phase (j) = tt (instead of 0) was 
performed. As expected, the ionization yields are a mirror 



image of the ones shown in Figs. [2] and [3j i.e. the lobes 
of the HOMO-2 point for <j) = tt toward positive z. The 
strong phase dependence found here even for a pulse as 
long as 8 cycles is very likely to be a consequence of the 
strong dipole moment of (oriented) water molecules. Thus, 
a very pronounced carrier envelope effect should occur for 
even shorter laser pulses. Since the dipole vector of the 
molecule is oriented along the z direction and the orbital 
density possesses mirror symmetry with respect to the yz 
and xz planes, the asymmetry occurs in the xz, but not 
in the xy plane. While the strong phase effect is interest- 
ing by itself and provides in sufficiently short pulses even 
information on the absolute orientation of a molecule and 
not only its orientation, it is problematic from the view- 
point of orbital imaging. In the case of the HOMO-2 of 
water it could appear as the correct image is obtained, but 
the "image" depends on pulse length and carrier envelope 
phase. Recall that in [1] the asymmetry was measured 
for HCl molecules using 50 fs pulses (and circular polar- 
ized light), but this required the coincident detection of 
ions in order to restore the molecular-frame information. 
Based on a simplified tunneling model, the found asymme- 
try was attributed to the dipole moment and not to the or- 
bital shape. Since a permanent or induced dipole moment 
influences the orientation-dependent ionization yield, this 
effect has to be corrected for, if orbital (densities) should 
be imaged. The size and importance of this correction in- 
creases with decreasing pulse length and thus, if a better 
time resolution should be obtained, by employing shorter 
pulses. 

The orientation-dependent ionization of the HOMO-1 
shows an asymmetry along the z direction similar to the 
one discussed for the HOMO-2. Correspondingly, the mir- 
ror image (with respect to the xy plane) is obtained for a 
phase (j) — t: and the asymmetry disappears for long pulses. 
Interestingly, a second asymmetry is observed when com- 
paring the cuts through the xz and yz planes, shown on 
the absolute scale in Fig. |4] The lobes in the xz plane are 
found to be substantially broader than the ones in the yz 
plane. Since the discussed difference in the widths is per- 
pendicular to the dipole moment (and shows the proper 
mirror symmetry) , it is clearly a consequence of the shape 
of the orbital. In agreement with the orbital shape, the 
ion yield in the xz plane is broader than in the yz plane, 
since the orbital density of the HOMO-1 is stretched in di- 
rection of the hydrogen atoms. However, for a sufficiently 
long pulse or an experiment with no control over the car- 
rier envelope phase inversion symmetry is automatically 
enforced and thus the widths of the lobes pointing in pos- 
itive or negative z direction will be equal in such a case. 
The different width is thus only a consequence of the upper 
lobe (pointing toward +z) of the orbital in Fig.|4] while the 
lower lobe {—z) dilutes the effect, since it is rotationally 
symmetric. 

This result is of interest, since it confirms directly that 
the orientation-dependent ionization refiects the position 
space orbital and not its momentum space counterpart, 
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Figure 4: (Color online) Shown are cuts of the ionization yield of 
the HOMO-1 in the xz and yz plane. The color coding agrees with 
the one in Fig. [2] Below, the shapes of the HOMO-1 projected on 
the xy, xz and yz plane are represented. 



since position and momentum spaces are reciprocal and a 
broad distribution in position space corresponds to a nar- 
row one in momentum space and vice versa. While this 
finding is in immediate agreement with the position-space 
based tunnel picture, it seems, on the first glance, to con- 
tradict the predictions of the strong-field approximation 
(SFA). Especially in its velocity-gauge version it predicts 
(within the frozen-core approximation) that the ionization 
rate of position-space orbital 4>j is related to its Fourier 
transform (kjv | 0j ) and thus to the momentum-space or- 
bital [Ml, 
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The momentum k^r of the emitted electron depends on the 
number of absorbed photons TV that is limited below by the 
minimum number A'o required to overcome the ionization 
threshold. The function / depends on the laser intensity 
/ and frequency w as well as on the projection of the mo- 
mentum of the emitted electron on the laser polarization 
vector e. For sufficiently high laser intensities, / peaks 
strongly at parallel emission of the electron [TT] and the 
integral over all emission directions reduces (very) ap- 
proximately to the integrand with kjv ~ e. Furthermore, 
if the laser frequency is low and the intensity sufficiently 
high, but not (much) higher than the over-the-barrier ion- 
ization threshold, the function / is strongly peaked at low 



momenta /cjv- Thus, according to SFA the low- momentum 
part of the electron density is probed. Low momenta corre- 
spond, however, to the electrons far away from the nuclei, 
ie. to the asymptotic part of the position-space orbital. 
Therefore, the predictions of the tunneling picture (prob- 
ing the electron density close to the barrier) and the ones 
from the SFA (probing the low-momentum component of 
the momentum-space orbital) are compatible and also in 
agreement with the present findings for H2O. Finally, in all 
formulations going beyond the frozen-core approximation, 
it is the Dyson orbital which should be more properly em- 
ployed. However, the Dyson orbitals of H2O are predicted 
to resemble the usual orbitals, at least for the HOMO, 
HOMO-1, and HOMO-2 considered in this work [39:. 

Figure |2] also shows the excitation probability of water 
for the given laser pulse. Again, the overall shape resem- 
bles the one of atomic p orbitals and is thus similar to the 
ionization yield. This similarity is most pronounced for 
the HOMO-1, while the HOMO excitation yield does not 
show the rotational symmetry around the x axis, but is 
compressed in the y direction. The HOMO-2 excitation 
yield shows on the other hand some additional structure 
that is absent in the ionization yield. The basic similar- 
ity of the orientational dependence of the excitation and 
ionization yields indicate the similarity of the processes 
leading to excitation and ionization as was discussed (for 
atomic helium) in [40j. A broad distribution of highly 
excited Rydberg states is populated, since (in the semi- 
classical picture) not all electron trajectories reach finally 
the continuum ( "frustrated tunnel ionization" ) . 

The structures in some of the excitation spectra and 
their absence in the ionization yields indicate on the other 
hand the occurrence of resonances that do, however, barely 
influence the ionization process. An exception is the 
HOMO-1 (small) structure in the yz plane that occurs 
close to the center. It is found in both the excitation 
and the ionization yields. Its origin becomes clearer from 
the comparison of the results that are obtained with the 
two different implementations of the SAE introduced in 
[T31 [3(2,. The results are shown in Fig. [5] Clearly, the 
structure occurs only in method B in which all orbitals are 
propagated in time, but is absent in method A in which 
all occupied field-free orbitals are excluded in the time 
propagation. In method A the ionization yields resemble 
in shape very well the one of atomic p orbitals. Note, 
also the small local minimum in the ionization yield of the 
HOMO that occurs along the y axis is absent in method 
A. The reason for these structures seen in method B may 
be interpreted as a coupling between the HOMO and the 
HOMO-1. As the HOMO is partially depleted due to ex- 
citation and ionization, it can partially be refilled from 
the HOMO-1. Due to its lower binding energy, this trans- 
ferred electron density is easier excited and ionized which 
leads to an enhanced excitation and ionization yield from 
the HOMO-1. At the same time, some population may 
be transferred from the HOMO into the HOMO-1 lead- 
ing to a stabilization and thus a reduced excitation and 
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Figure 5: (Color online) Comparison of the respective absolute 
main ionization yields of the HOMO Ifei (red) and the HOMO-1 3ai 
(blue) in the yz plane (top and middle panels) and in the xz plane 
(bottom panel). Laser parameters are the same as in Fig. [5] 0° mean 
alignment of the electric field polarization along the z axis and 90° 
along the x or y axis. The two upper graphs differ by the method 
applied: Method A (top) and method B (middle). 



Figure 6: (Color online) Ionization yields in the xz (left) and yz 
plane (right) for the HOMO and HOMO-1 for a 800 nm, 20 cycles 
pulse at the intensity 5 X 10^^ W/cm^. Our SAE-TDSE results are 
compared to TDVFD results [21] • The TDVFD data is scaled (by 
the factors given) to agree with the SAE-TDSE data at 0°. 



ionization from the HOMO. A similar coupling of initially 
occupied (field- free) orbitals was found for CO2 T^, but 
there the influence on the ionization yield was much more 
pronounced. It should be understood, however, that this 
interpretation is given in terms of field- free eigenstates that 
are simply basis states in the TDSE calculation, and thus 
loose a direct physical meaning during the pulse. 

For the visibility of the different orbitals in an experi- 
ment the absolute magnitude of the different ion yields is, 
of course, important. As is evident from Fig. [5] within the 
xz plane the HOMO-1 ionization yield is much larger than 
the one from the HOMO due to the nodal plane of the lat- 
ter. On the other hand, the space volume in which the 
HOMO-1 ionization yield dominates, is rather small, since 
it extends only to about 10% into the yz plane. For the 
laser pulse adopted here, a very well aligned H2O sam- 
ple would be required in order to observe HOMO-1 ion- 
ization. Integrated over all directions ionization is domi- 
nated by the one from the HOMO (note the logarithmic 
scale). A dominant ionization from the HOMO-1 in the 
xz plane had been predicted in [2T] based on a TDVFD 
calculation. A direct comparison (with identical 20 cycle 
pulses with peak intensity 5 x 10"'^^ W/cm^) is shown in 
Fig. [6] Qualitatively, the results are very similar, though 
the present ones tend to be broader, i. e. minima and max- 
ima are slightly less pronounced. On the absolute scale, 
the results differ by factors in between 1.7 and 10. 

It is also interesting to consider the intensity depen- 
dence of the relative ionization yields of the HOMO and 
the HOMO-1 shown in Fig. [7) For low intensities, the max- 
imum ionization occurs along the direction of the HOMO 
(y axis), followed by the one of the HOMO-1 along its 
maximum extension (z axis) . The ionization of the HOMO 
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Figure 7: (Color online) Ionization yields of ttie HOMO Ifei (red) 
and the HOMO-1 3ai (blue) for alignments of tlie electric field vector 
E along the three main axes x (squares), y (circles) and z (diamonds) 
for 8-cycle cos-^ laser pulses (800 nm). All ionization yields from 
the SAE-TDSE calculation were multiplied by a factor 2. The inset 
shows the ratio of the HOMO ionized along the j/-direction (red filled 
circles) and the HOMO-1 ionized along the z-direction (blue filled 
diamonds) together with an ADK estimate of the ratio. The arrows 
indicate the intensity used in most other graphs, especially in Fig. [2] 



along the x and z directions is practically identical and, 
for not too high intensities, agrees also to the one of the 
HOMO-1 along the y direction. The lowest probability is 
found for the HOMO-1 along the x direction, due to the 
absence of the resonant structure observed in y direction 
and discussed above. Comparing the ionization probabili- 
ties of the HOMO and the HOMO-1 along the directions 
of maximum ionization (inset of Fig. [?]) one finds a pro- 
nounced intensity dependence. Clearly, at high intensities 
the ratio approaches 1, since the HOMO ionization satu- 
rates earlier than the one of the HOMO-1. In a comparison 
to experiment, saturation requires to consider also focal 
volume effects. Furthermore, the transition from a factor 
2 to 1 due to saturation of the single-electron ionization 
should differ for HOMO and HOMO-1, if it follows the 
earlier findings for H2 where a rather sharp transition be- 
tween factors 1 and 2 for about 10 % ionization probability 
was found. 

In a recent experiment the high-harmonic spectra of 
H2O and D2O have been measured. In the analysis based 
on a calculation using the present theoretical approach it 
was concluded that the harmonics for the different isotopes 
differ significantly, because ionization from the HOMO-1 
excites a vibrational wavepacket, while HOMO ionization 
does not [5] . In view of the present finding it is clear that 
the observability of the HOMO-1 high harmonics in an 
isotropic sample requires that the lower ionization proba- 
bility is partially compensated for by the higher recombi- 
nation rate (3rd step of the three-step model). Further- 
more, as is also seen from the present results, saturation 
plays an important role for higher intensities, since the ra- 
tio of the maxima of the ionization yields from the HOMO 
and the HOMO-1 approaches 1 for high intensities (cf. the 
inset of Fig. [?]). An estimate of this ratio with ADK (in- 
cluding the saturation effect) is also shown. Clearly, ADK 
overestimates the ratio significantly for lower intensities 
and reaches the ratio 1 due to saturation very abruptly. 
Concerning the direct imaging scheme discussed above, the 
intensity dependence of the ratio is important, since the 
image of the HOMO density may be distorted by an ad- 
mixture of lower lying orbitals and this admixture changes 
with intensity. While this is very critical for a measure- 
ment at a single intensity, the intensity dependence of the 
ratio can, on the other hand, also help to disentangle the 
contributions from different orbitals. 



6. Conclusion 

In conclusion, the response of water molecules to ultra- 
short intense laser pulses (at 800 nm) was investigated by 
solving the time-dependent Schrodinger equation within a 
single-active-electron approximation, but taking into ac- 
count a complete multicenter molecular core. In fact, the 
influence of the laser field on all electrons was consid- 
ered, but the electron-electron interaction is restricted to a 
single-determinant approximation. It was found that the 
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orientation-dependent ionization from the three energeti- 
cally highest occupied orbitals reflects their overall shape. 

However, for a sufficiently long pulse the image acquires in- 
version symmetry, since the electric field vector points ap- 
proximately symmetrically up and down. For short pulses, 
in the present case even as long as 8 cycles, deviations from 
the inversion symmetry of the ionization distribution were 
found. Although it may be a consequence of the asymmet- 
ric shape of the orbital, the visibility of this effect even for 
such relatively long pulses suggests that it is mainly due 
to the permanent dipole moment of an (oriented) water 
molecule. While this effect provides additional informa- 
tion, it distorts the image of the orbitals and it depends 
evidently on the absolute phase and time duration of the 
laser pulse. On the other hand, the missing rotational 
symmetry of water compared to the previously mainly dis- 
cussed linear molecules allowed a distinction of broad and 
narrow fcatTircs of the orbitals, provided they survive the 
enforced inversion symmetry. 
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